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Abstract 



£>"v Eternally inflating universes can contain large thermalized regions with dif- 

ferent values of the constants of Nature and with different density fluctuation 
spectra. To find the probability for a 'typical' observer to detect a certain 
set of constants, or a certain fluctuation spectrum, one needs to compare the 
volumes occupied by different types of regions. If the volumes are taken on 



CO 

o 

■ an equal-time hypersurface, the results of such a comparison are extremely 

o 
in 

On 



sensitive to the choice of the time variable t. Here, I propose a method of 



comparing the volumes which is rather insensitive to the choice of t. The 

i 

^jy method is then applied to evaluate the relative probability of different min- 

ima of the inflaton potential and the probability distribution for the density 
^ ■ Gucluat ion sped ra. 



I. INTRODUCTION 

Inflation is a state of rapid (quasi-exponential) expansion of the universe [l]. The infla- 
tionary expansion is driven by the potential energy of a scalar field (p (called the 'inflaton'), 
while the field slowly 'rolls down' its potential V(<p). When tp reaches the minimum of the 
potential, this vacuum energy thermalizes, and inflation is followed by the usual radiation- 
dominated expansion. 

One of the striking aspects of the inflationary cosmology is that, generically, inflation 
never ends @|§. The evolution of the field <p is influenced by quantum fluctuations, and 



as a result thermalization does not occur simultaneously in different parts of the universe. 
Inflating regions constantly undergo thermalization, but the exponential expansion of the 
remaining regions more than compensates for the loss, so that the inflating volume keeps 
growing with time. 

Eternally inflating universes may contain large thermalized domains with different values 
of the parameters determining the low-energy physics. For example, if gravity is described 
by a Brans-Dicke-type theory, different domains may have different values of the effective 
gravitational constant Q. Another possibility is that the potential V(ip) has several minima 
with different masses and couplings of light particles (and different symmetry breaking 
schemes). Although we cannot observe other regions with different low-energy physics, a 
study of their properties and their relative abundances in the Universe is not necessarily 
pointless. If we assume that we are a 'typical' civilization inhabiting the Universe, then such 
a study may help us understand why we observe the particular set of low-energy constants in 
our part of the universe |§-§]. (The assumption of being typical was called the 'principle of 
mediocrity' in Ref. |J. It is a version of the 'anthropic principle' which has been extensively 
discussed in the literature ||.) 

A related question was raised in a very interesting recent paper by A.Linde, D.Linde 
and A.Mezhlumian |7j who studied the spectrum of density fluctuations seen by a 'typical' 
observer in an eternally inflating Universe. They came to a surprising conclusion that such 
an observer should find herself near the center of a very deep minimum of the density field. 
If correct, this result may rule out a wide class of inflationary models. 

One could try to implement the mediocrity principle by comparing the number of civ- 
ilizations inhabiting regions with different values of the constants (or different fluctuation 
spectra) at a given moment of time. This was the approach taken by Linde et.al. 0,|6]- 0. 
However, one finds that the resulting probability distributions are extremely sensitive to the 
choice of the time coordinate t [§,|J. Since no particular choice of t appers to be preferred, 
this casts serious doubt on any conclusions reached using this approach. In the present paper 
I attempt to define the probability distribution for the constants in a coordinate- independent 
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way. 

The idea is to think in terms of the total number of civilizations in the entire spacetime, 
rather than their number on a particular spacelike hypersurface, t = const. To be specific, 
let us consider 'new' inflation with a potential V((p) of the form illustrated in Fig.l. The 
potential has two minima, and the values tp^ and (p^ near the minima correspond to the 
end of inflation. The relative probability for a civilization to be in one type of thermalized 
region versus the other can be estimated as |J 

V (2) ~ V (2) V (2)" ' ' 

K * "civ 



Here, V* is the 3-volume of the hypersurface(s) <p = <pi\ j = 1,2, and u^l is the average 
number of civilizations originating per unit volume (I assume that civilizations can 

originate for only a finite period of time after thermalization, so that their number per unit 
thermalized volume is finite). The dependence of the 'human factor' v c i v on the low-energy 
physics constants has been discussed || in relation to the anthropic principle. Here, I will 
disregard it and concentrate on determining the ratio 

r = V«/V< 2 >. (2) 



In models where the Universe is closed and inflation is not eternal, the volumes V: 



U) 



are finite, and the ratio r is well defined. However, in an eternally inflating universe, 

(2) 

and V* are infinite and have to be regulated. If one simply cuts them off by introducing 
a hypersurface t = const and including only parts of the volumes in the past of that hy- 
persurface, then again one finds that the ratio r is highly sensitive to the choice of the 
cutoff hypersurface, so we are back to our original problem. Here, I would like to suggest 
an alternative regularization procedure which appears to give more reasonable results and 
which is more robust with respect to changes of the time variable. 

The rest of the paper is organized as follows. The general spacetime structure of an 
eternally inflating universe is discussed in the following Section. The proposed regularization 
procedure is outlined in Section 3, and a formalism for calculating relative probabilities for 



different minima of V(ip), which is based on this procedure, is introduced in Section 4. 
The dependence of the resulting probabilities on the choice of the time variable is discussed 
in Section 5, with a specific example worked out in Section 6. The spectrum of density 
fluctuations detected by a typical observer is calculated in Section 7, and the result is 
compared to that of Linde et.al. |7J]. The main results of the paper are summarized and 
discussed in Section 8. 

II. INFLATIONARY SPACETIMES 

An inflating Universe can be locally described using the synchronous coordinates, 

d s 2 = dr 2 -a 2 (x,r)dx 2 . (3) 

The lines of x = const in this metric are timelike geodesies corresponding to the worldlines 
of co-moving observers, and the coordinate system is well defined as long as the geodesies 
do not cross. This will start happening only after thermalization, when matter in some 
regions will start collapsing as a result of gravitational instability. Hence, the synchronous 
coordinates @ can be extended to the future well into the thermalized region. 

The classical evolution of the scale factor a(x, r) and of the scalar field y(x, r) is described 
by the equations 

(a/a) 2 = H 2 ^ 8ttV(^)/3, (4) 

ip « -H'(tp)/4ar, (5) 

where dots represent derivatives with respect to r and I use Planck units, % = c = G = 1. 
The potential V(<p>) is assumed to be of the form shown in Fig.l. It is assumed also that 
the field <p is a slowly-varying function of the coordinates, so that spatial gradients of ip 
can be neglected and ip 2 <C 2V(cp). With the aid of Eqs.(|4]),(|), the latter condition can be 
expressed as 
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H' < 6H. (6) 

This condition is violated near the points p> = <p*\ signalling the end of inflation. The slow 
variation of (p implies that H is also a slowly-varying function of x and r, and thus the 
spacetime is locally close to DeSitter, with a horizon length if -1 . 

Quantum fluctuations of the field ip can be pictured as a 'random walk' [superimposed 
on the classical motion (gj)] in which tp undergoes random steps of rms magnitude (5(p) rms = 
H/2tt per Hubble time, 5t = H , independently in each horizon-size region (£ ~ H -1 ) ||10|| . 
The fluctuations are dynamically unimportant if the classical 'velocity' \(p\ is much greater 
than the characteristic speed of the random walk, (Sip) Tms /5r = H 2 /2ir, which gives 

H'^>H 2 . (7) 

This condition is violated in the region c^q 1 ^ < tp < tpQ near the top of the potential 
(see Fig.l). The dynamics of tp in this region is dominated by quantum fluctuations. The 
deterministic slow-roll regions are bounded by the values tp^ and tp* . 

The spacetime of the inflating universe is schematically represented in Fig.2. Thermal- 
ized regions of different type are shown by different kinds of shading; their boundaries are the 
surfaces tp = tp* . These boundaries expand into the inflating region, but become asymptot- 
ically static in the coordinate space as r — ► oo. (If the speed of expansion remained finite at 
large r, then thermalized regions would merge and inflation would terminate everywhere). 
This does not mean, however, that the boundary surfaces become timelike at large r. To 
picture the geometry of these surfaces, one has to keep in mind that physical lengths in 
an inflating universe are related to coordinate differences by an exponentially growing scale 
factor. As a result, the boundaries of thermalized regions (as well as all surfaces of constant 
tp in the slow-roll regime) are spacelike, and are in fact very flat. 



The spacelike character of the constant-y? surfaces can be understood as follows ||TT 
Consider the normal vector to the surfaces, d^ip. We have 

d^ip = tp 2 - a- 2 (Vtp) 2 . (8) 



The spatial gradients of p are caused by quantum fluctuations. On the scale of the horizon, 
the gradient is of the order a -1 |V<£>| ~ (Sp) Tms / H' 1 ~ H 2 /2ir, and is even smaller on larger 
scales. On the other hand, from Eqs.(^),(^), \p\ ^> H 2 /2tt. Hence, d^pd^p > 0, and 
the surfaces <p = const are spacelike. Moreover, these surfaces are locally nearly parallel 
to the surfaces r = const (which correspond to horizontal lines in Fig.2). This can be 
seen by considering the scalar product of the unit normals, = (1, 0, 0, 0) and = 

d^id^d^)- 1 ' 2 , 

w' i Ti A4 = [l-o- 2 (V^)7^]- 1 / 2 «l. (9) 

For an internal observer in one of the thermalized regions, the surface <p = ip* at the 
boundary of that region plays the role of the big bang. The natural choice of the time 
coordinate in the vicinity of that surface is t = <p, so that the constant-t surfaces are 
(nearly) surfaces of constant energy. Since these surfaces are infinite, the observer finds 
herself in an infinite thermalized universe, which is causally disconnected from the other 
thermalized regions. The situation here is similar to that in the 'open-universe' inflation 
|T2|] where thermalized regions are located in the interiors of expanding bubbles and have 
the geometry of open (k = —1) Robert son- Walker universes. 

We see that, depending on one's choice of the time coordinate t, an equal-time surface 
may cross many thermalized regions of different types (e.g., for t — r), may cross no such 
regions at all (say, for t = p with <p in the slow-roll range), or may lie entirely in a single 
thermalized region. Not surprisingly, with a suitable choice of t, one can get any result for 
the volume ratio r in Eq.(0). 

III. THE PROPOSAL 

Let us choose a region in which the field <p is close to the top of the potential (<p ~ 0) at 
some moment of time in the coordinates ([$]). We can set r = and a = 1 at that moment. 
The choice of the initial moment and of the particular region is unimportant, due to the 
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self-similar nature of eternal inflation Let us now introduce a large number N of 

co-moving observers uniformly spread over the region. Some of these observers will end up 
in thermalized regions of type 1 and others in regions of type 2 (see Fig. 3). The initial 
3-volume of the region Vo has to be sufficiently large, so that the corresponding co-moving 
volume includes a large number of thermalized regions of both types. (Alternatively, one 
could consider an ensemble consisting of a large number of small regions). The initial volume 
per observer is Vq/N, and if a*j is the scale factor for the z-th observer at the time when she 
reaches </?* , then the 3- volumes of the hypersurfaces tp = <pi can be represented as 

Cf) 

Vi^W 1 ^, (10) 

i 

Here, j = 1, 2, and the summation is over all observers ending up in thermalized regions of 
type j. [I disregard the factors (w^n^) -1 ~ 1; see Eq.@]. 

The sums in Eq.([l0|) are, of course, divergent in the limit N — > oo. My proposal is 
to regularize these sums by cutting off a small fraction e of terms having the largest scale 
factors a*, with the same value of e for both types of thermalized regions. The volume ratio 
r in Eq.(0) can then be evaluated in the limit e — > 0. 

A similar regularization can be implemented using the proper time r, or some other 
'time' t, as a cutoff variable instead of the scale factor (that is, discarding a fraction e of 
observers having the largest values of t when they reach (p = (p*^). We shall see that the 
volume ratio r is not very sensitive to the choice of the cutoff variable. 

The reader may be concerned that coordinate-dependence in this approach is introduced 
from the very beginning, when we choose the initial volume Vo on an equal-time surface 
in the synchronous coordinate system. The role of this volume, however, is only to define 
an appropriate congruence of geodesic observers. Geometrically, the construction can be 
described as follows, (i) Find a spacetime region where <p ~ 0. (ii) Choose a hypersurface 
in this region which has approximately constant, small intrinsic and extrinsic curvature 
(^R « const ;$ H 2 ). (iii) Consider a congruence of geodesic observers whose worldlines are 
orthogonal to this hypersurface. 



IV. CALCULATION OF THE VOLUME RATIO 



With the above prescription, the volume ratio r can be calculated in terms of the prob- 
ability distribution V(p, t) for the inflaton field p. This distribution satisfies the diffusion 
equation PJHI 



d t V = -d v J, 



(11) 



where the flux J(p>, t) is given by 



(12) 



The first, 'diffusion' term in Eq.([TJ) represents quantum fluctuations of the field p, while the 
second term corresponds to the classical 'drift' due to the potential V(ip). The parameter 
a is defined so that for a = 1 the time variable t is the proper time r, while in general t is 
related to r through 



dt = H x - a dT 



(13) 



In particular, for a — 0, t — In a. The boundary conditions at p = (p* are 



J(p* 



(47r)- L H a - L H'V(ip*). 



(14) 



The initial distribution V(<p, 0) is sharply peaked near p = and is normalized so that 



S 2 > 



V(p,0)dp = l. 



(15) 



The implementation of the e-cutoff procedure is most straightforward with the choice of 
a = 0. Then a — e*, and Eq. (|10|) gives 



V? J = V 



J(pi^\t)e 3t dt 



(16) 



where the cutoff 'times' ti^ are determined from 



U) 



J(<p<?\t)dt 



J(<p®,t)dt 



ep 



U) 



(17) 



The conservation of probability implies that 

p(i) + p (2) = ! (18) 

Equations (|T5|) , QT 7 ?!) are easily understood if we note that | J(<p*\ t)\dt is the fraction of all 
observers who end inflation in a region of type j at time t in the interval dt. We note also 
that the cutoff times — > oo as e — ► 0. 

The solution of Eq. flTTD can be represented in the form 

oo 

V{^t) = Y.M^)e- ln \ (19) 

n=l 

where < 71 < 72 < .... In the limit t — ► 00, 

7%, t -> 00) = ^i(^)e" 7lt . (20) 

The eigenfunction ^1 (<£>) can be found explicitely in the slow-roll regions, away from the 
top of the potential, where the diffusion term in Eq.(|T2D is negligible. The equation for 
i/ji(<p) in these regions takes the form 

d^H^H'^) = -47T 7l Vi, (21) 

and its solution is 

... TJ(m\ T Ff> 1 

(22) 



The solution ipi\ip) applies in the range between (pffl and (p*\ [Note that Eq.(p2[) satisfies 
the boundary condition (0)]. The coefficients c^- 1 and the eigenvalue 71 can be found by 
matching Eq.(|22|) to the solution for ipi in the range ip^ < (f < <f^. An example will be 
given in Section 6. 

We shall now use the asymptotic forms (f20|) , (f22"D to evaluate the volume ratio r. For 
inflation to be eternal, the exponential decay of the distribution function (^) should be 
slower than the expansion of the inflating regions, a 3 = e 3t , so that the total inflating 
volume grows. Then, 71 < 3, and the integral in Eq.([R]) for V* is dominated by the upper 



limit. With the aid of Eqs. (|20|) , (]22|) we can write 



yO') = [4vr(3 - 7i)] _1 V c (j) exp [(3 - lx )t® + 7l t<-» 



(23) 



where 



An / — ^-df 



(24) 



is the 'time' of the classical slow roll from tp^' to (p* [see Eqs.(|5|), (|l3l) 1. 

The cutoff 'times' t^p are determined by Eq.([l7|), and substituting the asymptotics 
(13) '(H) i n ^ ne left-hand side of that equation, we have 



(4tt 7i )- 1 c {j ' ) exp[ 7l (t^ - tf)] = ep U) . 



(i) 



(25) 



Finally, substituting this into (|23|), we obtain 

V (j) = c (j)v [4tt(3 - 7l )]- 1 [4vr 7l ep (j) /c (i) ]~ (3 ~ 7l)/71 exp[3^] 



(26) 



and 



(i) 



V. 



(2) 



~(1) 



3 /7i r^(2)"| ( 3 ~7l)/7i 



(1) 



Here, 



exp 



P 



-4vr / ,.. ^±d£ 



i* 

7M 



(27) 



0) 



(28) 



is the expansion factor during the classical slow roll from tp^ to <p* . 

The quantities 7 i, p^\ and Z* all depend on the shape of the potential V(^). The 
eigenvalue 7 i and the ratio p^ /p^ are determined mainly by the shape of the potential near 
p = 0, where the dynamics of tp is dominated by quantum fluctuations. Unless V(ip) is very 
asymmetric near ip = 0, the numbers of observers ending up in the two types of thermalized 
regions will have the same order of magnitude, and thus pW/pW ~ 1. If the potential is 
approximately symmetric near the top, V(— p) ~ V(p>), then we can set (p^' « —p>^\ and 
one expects that ■?/>i(y 3 o 1 ' ) ) ~ V ; i(v 9 o 2 ' ) )' that is, c^/c^ ~ 1. In this case, the most sensitive 
dependence of r on the shape of V(ip) is through the expansion factors Z* , 

r-fzW/Zf 2 )] 3 . (29) 
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V. DIFFERENT TIME VARIABLES 



With a different choice of the time variable (a 7^ 0), the calculation of r is somewhat 
more complicated. The difficulty is that the 'time' t is no longer related to the scale factor 
a, and we cannot write the analogue of of Eq.flTED for the volume V* . This can be done 
by introducing another distribution function V(<p,t), such that V(ip,t)dip is the physical 
volume occupied by regions with ip in the interval dip at time t. In other words, V(ip,t)dip 
is the 3-volume on a hypersurface t = const in which the inflaton field has values between 
ip and ip + dip. The function V((p, t) satisfies the modified diffusion equation [15|-[17|,|| 



d t V = -cLJ + 3H a V, 



(30) 



where J is given by Eq.(|T2"D with V replaced by V. At the initial moment, t — 0, the 
distributions V and V are proportional to one another, 



V(<p, 0) = V V^,0). 



(31) 



The volume in which <p> rolls down to the value <p£ in a time interval dt is \J{<p ( * 3> , t)\dt, 



,0") 



and thus the volume V*^ with a cutoff at can be written as 



V (i) 



J(<pP,t)dt 



(32) 



The cutoff times ^ are still determined by Eq. (|T7D . Hence, in order to calculate r for a 7^ 0, 
one has to solve both equations ( [TT]) and ( |30| ) for P and P. 
The distribution V((p,t) can be expanded in eigenfunctions, 



(33) 



n=l 



where 71 > 72 > •••• For inflation to be eternal, we should have 71 > 0. In the limit t — > 00, 



v{<p,t^ 00) = M^ 11 . 



(34) 



As in the previous Section, the form of ipiif) i n the slow- roll regime can be found by 
neglecting the diffusion term in Eq.(|30|). 



11 



d^H^H'^x) + 127rH a iPx = 4tt7i^i. 



(35) 



The solutions of this equation in the regions between yjj, and ^7* 7 '' are 



cxp 



f¥> #(£) /•¥> H 1 - ^) 



(36) 



with c^) = const, and combining Eqs.(|32]),(|3"4"D and (|36|), we have 
Here, Z* is given by Eq . fl28|) and has the same meaning as in Eq. 



t[ j) = -47T / 



U) 



(37) 

and is given by 
(38) 



[Note that for a = 0, Eq.(|35|) reduces to Eq.(pl|) with 71 = 3 — 71. Then it is easily 
understood that ffl = V c {j) , and thus Eqs.(H),(|D agree with Eqs.©,(|7|)]. Now, with 
the cutoff times from Eq . (p5|) , we obtain 



Vi 



(1) 



(2) 



c (l) p (2)] 7l/7i 



'* 

R2j 



(39) 



c (2)p(l) 

The coefficients c^, and the eigenvalues 71, 71 depend on the choice of the parameter 
a, and thus the ratio r has some dependence on the time parametrization. [The quantity 
is the fraction of the observers that end up in the j-th type of thermalized regions, and 
should be independent of a]. However, this dependence appears not to be very strong. In 
particular, in cases where the potential V(ip) is nearly symmetric in the region dominated 
by quantum fluctuations near the top, we have c^ l > / ~ /c^ 2 ' ~ pi 1 ) /pi 2 ) ~ i ; and the 
estimate (|29|) holds for any value of a. 

In contrast, the result would be extremely sensitive to the choice of a if we used a cutoff 
at a fixed value of t. Then we would get 





5(1) 


>(!)" 




~ C( 2 ) 


7(2) 



expHy^M-tf)] 



(40) 



The quantity is the time it takes the scalar field (p to roll from </?q''' to tp*' . It obviously 
depends on the choice of the time coordinate. 
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0) 



,0") 



To assess the importance of the last factor in (|40|), we need an estimate for 71. It can be 



shown that 18 



7i < 3flSL», (41) 

where H max is the largest value of H{tp). [H max = H(0) for the potential in Fig.l]. For a 
generic potential, 

71 ~ H^ ax (42) 
gives a reasonably good order-of-magnitude estimate. The quantity 

d = ^H- a ax < 3 (43) 
can be called the fractal dimension of the inflating universe [I3,iT5|. The values of d <C 1 can 



be obtained only when the potential is fine-tuned so that eternal inflation is barely possible. 
Returning now to Eq.(|40|), we see from Eqs.([28l),(|3~8"l), and (f42[) that, if a > 0, then 



7itl J ^ > 31n.z¥\ and thus the last factor in ( |40D is typically no less important than the ratio 
of the expansion factors. In particular, cutoffs at a fixed proper time (a = 1) and at a fixed 
scale factor (a = 0) often give drastically different results. 



VI. AN EXAMPLE 

The eigenfunctions ipi((f), ^(v 9 ) (and the corresponding eigenvalues) can be approxi- 
mately found in the whole range of (p for a potential V(<p) shown in Fig. 4. The poten- 
tial is completely flat (V = Vq) for (p^ < <p < tp^\ satisfies the slow-roll conditions 
outside this range, and is arbitrary otherwise. Without loss of generality, we can 

+. ( 2 ) (1) — 

set ip'o = -ip K = <p . 

In the flat region of the potential, \<p\ < <p , the equation for V ; i(v 9 ) is 

4>'( + k 2 ^ = 0, (44) 

where 
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K 2 = 87r 2 V a+2) W-7i) ( 45 ) 
and H = (87TV0/3) 1 / 2 is the expansion rate in the flat region. The solution of Eq. fl44|) is 

ipi(<f) = A\ cos(k<^) + A 2 sm(K<p). (46) 



In the slow-roll regions at \cp\ > <po we can use the approximate solutions (p6|). The match- 
ing onditions at if = ±<po require that ipi((p) and the corresponding flux J\{(p) should be 
continuous, 

[^Ww, = [Jl\<P=±<Po = °- ( 47 ) 

Combining Eqs. d36|),(^6|) and (0), we obtain the following system of equations 

A l cos(^o) - A 2 sin(^o) = c (1) H^ a / H[, (48) 

A x cos(/?v9 ) + A 2 sin(^o) = c {2) H^ a / H' 2 , (49) 

kH^ 2 ^! sin(^ ) + A 2 cos(k^o)] = 2tt5 (1) , (50) 

kH^ +2 [-A 1 sin(^o) + A 2 cos(«p )] = 2tt£ (2) , (51) 
which yields, after some algebra, 

Pi- p 2 + 2K(fi = nir, (52) 



^. ta (A±|±=). ,53, 



^ , n cos/3 ; 



cos (3 2 

Here, I have introduced the notation 



(54) 



Pj = tan' 1 (uHl/ZiiH'j), (55) 
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where Hj is the derivative H'(<p) taken on the slow-roll side of the point ip^ . 

Since H[ > and H' 2 < 0, the left-hand side of Eq.fl52|) is a monotonically growing 
function of k. This equation, therefore, has a single solution for any value of n. From 
Eq.(f^5"D it is easily seen that the largest eigenvalue 71 corresponds to the smallest value of 
k 2 , and thus to n = 1. 

For eternal inflation to be possible, we need 71 > 0, and Eq.(IS) yields k 2 > 24tt 2 Hq 2 . 



Together with the condition (|7J) this implies \/3j\ <C 1, and from Eq.(^) with n = 1 we have 

K 7H 7r/2(/9 . (56) 

The condition for eternal inflation to be possible in this model can now be written from 



Eqs.(H) 



> H Q /AVQ. (57) 



The meaning of this condition is that the size of the flat region 2(po should not be smaller 
than the size of the random walk 'step' Hq/2tt. 

The equations for the function V ; i(v 9 ) are essentially identical to those for V ; i(v 9 )- Apart 
from the trivial replacements -^l - ► ifti, — > the only change is that Eq.(f45l) is replaced 
by 

k 2 = 87r 2 F - (Q+2) 7i. (58) 



Since the smallest eigenvalue 71 again corresponds to the smallest solution of Eq.flSl^) for k, 
we conclude that the value of k is the same for both functions, 4>i((p) and 4>i(<p). Hence, 

71 24tt 2 



71 H 2 k 2 

and 



1 (59) 



c (2) /c (l) = g (2) /g (l)_ (6Q) 

Now we are ready to analyze the dependence of the volume ratio r in Eq.(^) on the 
parameter a. This parameter does not appear in Eq . (|52|) for k, and thus k is independent 
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of a. Then it follows immediately from (|59|), ( |60D and ( ^4|) that 71/71, & 2 > /o> x \ and c^'/d 1 ' 
are also independent of a. Since and are generally a-independent, we conclude that, 
for the model considered in this Section, the ratio r is completely insensitive to the choice 
of time parametrization. 

This conclusion may seem somewhat surprising, but in fact it is not difficult to under- 
stand. A special property of our model is that there is a one-to-one correspondence between 
the time t when a co-moving observer reaches (p = (p* and the expansion factor a along her 
worldline at that time, 

a = exp[HZ(t-t<p)]zP. (61) 

For this reason, we can use a time coordinate with any value of a as a cutoff variable, and 
the result will be the same as with the scale factor cutoff. 



VII. DENSITY FLUCTUATIONS 



Let us now consider the spectrum of density perturbations detected by a 'typical' ob- 
server. These perturbations are determined by quantum fluctuations of the inflaton field <p 



along the observer's worldline [^UJ. For an ensemble of observers at the same value of H, 
the fluctuations 5(p are random Gaussian variables with a distribution 



dV (S V ) = (27Rx)- 1/2 exp 



2<X 2 



(62) 



where 



a = H/2 



71 . 



(63) 



On a given co-moving scale, the fluctuation is produced (that is, it becomes a part of the 
classical field) at the time when that length crosses the horizon, and the value of H in 
Eq.(|63"D should be taken at that moment. (For the purposes of this discussion, it is sufficient 
to take a simple-minded view that the fluctuation occurs instantly). The gauge-invariant 
density fluctuation on the corresponding scale is 



16 



Sp/p ~ 87rH5(p/H', (64) 
where H' = dH/dip and again all quantities are taken at horizon crossing. Averaging over 



the distribution (p2n gives the standard result [20 



{Sp/p) rms ~AH 2 /\H'\. (65) 

Linde, Linde and Mezhlumian [|7| have pointed out that the distribution (|62"D is not 
necessarily identical to the probability distribution dV(5<p) for the values of Sip detected by 
a 'typical' observer. There may be a correlation between Sip and the expansion factor a, 
in which case observers with different values of Sip will contribute to dV(Sip) with different 
weights. Linde et.al. arrived at a surprising conclusion that 'typical' fluctuations of ip can be 
much greater than the rms value (Q). Here, I shall first calculate dV(5ip) using the e-cutoff 
procedure and then compare it with the results of Linde et. al. 

Fluctuations of <p on different length scales are statistically independent and can be 
treated separately. We can therefore concentrate on a single fluctuation at some value 
ip = ifi (the same for all observers), disregarding all the rest. We shall assume that this 
value is in the slow-roll regime. 

To find the probability distribution for Sip, we start with the ensemble of co-moving 
observers defined in Section 3. We can split it into sub-ensembles, such that the fluctuation 
Stp at ip — ip\ is the same with an accuracy dSip for all observers in each sub-ensemble. The 
numbers of observers in different sub-ensembles are proportional to dVo(Sip). The probability 
dV(Sip) is then given by 

rfP(M«E4, (66) 

i 

where the summation is done over the observers in the corresponding sub-ensemble. (Here 
we consider only observers ending up in the same type of thermalized regions). The sum 
is regularized by discarding a fraction e of observers with the largest values of a*. (A more 



careful analysis shows GBJ that the same results are obtained if the sum is cut off using a 



time variable with an arbitrary value of a). 
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When the field p> undergoes a quantum jump Sip, the amount of proper time r necessary 
to complete the classical rollover to is changed by 

St = Air5(p/H[, (67) 

where the subscript '1' indicates that the corresponding quantity is taken at <p — <pi. As a 
result, the volume factor acquires an additional factor 

f(5ip) = exp(3H 1 Sr) = exp[12ir(Hi/H[)5(p]. (68) 



Since this factor is the same for all observers of the sub-ensemble, the cutoff in the sum 
is not affected by the fluctuation. That is, the discarded observers are the same as would 
be discarded for Sp = 0. 

Apart from the overall number of observers and the magnitude of the fluctuation at ip = 
<px, different sub-ensembles are statistically equivalent. Hence, the probability distribution 
dP(5<p) is proportional to 



dV(Sp>) oc dVo(Sp)f(Sp) oc exp 



d(5<p). (69) 



The distribution (|6~9D describes fluctuations with a non-zero mean value, 

< Sip >= 3Hf/nH[. (70) 

Deviations from this mean value are still Gaussian, with the dispersion 

a = [< {Sip- < Sip >f >] 1/2 = HJ21T. (71) 

According to Eq . (|70|) , the inflaton tends to fluctuate in the direction opposite to the classical 
roll [see Eq.([|)]. This is easy to understand: 'backward' fluctuations prolong inflation and 
increase the expansion factor. From Eqs.(|64D,(|70D, the mean density fluctuation is 

< Spi > jp ~ 24Hf/H[ 2 . (72) 

The average density profile < 5p(x) > / p can be easily found for a given H(ip). 



The appearance of a non-trivial average density distribution is an interesting, and in 
principle observable, effect. However, the magnitude of the effect is hopelessly small. From 
Eqs.(|65|) and ( |72|) we have 



On scales of astrophysical interest, (Sp/p) rms ~ 10~ 5 , and (assuming the observed density 
fluctuations are due to inflation) < 5p > / p ~ 1CT 10 . 

These results are very different from those obtained by Linde et.al. who studied the 
probability distribution for 5(p on equal-time surfaces, t = const. Using the proper time r 
as the time coordinate, they considered an ensemble of observers who reach (p = (p* at a 
given moment, with the weight assigned to each observer. They concluded that a typical 
observer will detect large quantum jumps of <p> in the direction of the classical roll. This 
conclusion can be qualitatively understood as follows. If I have to reach the point <p = ip* 
at a given time r with the largest possible scale factor a*, then the winning strategy for me 
is to spend as much time as possible at values of (p where V(<p) is large and the expansion 
is fast, and then quickly rush towards (p*. 

To make this quantitative, it is convenient to consider fluctuations that bring (p to a given 
value (p\. Any further fluctuations can be disregarded due to statistical independence. The 
evolution between <p\ and is the same for all observers. In particular, it takes the same 
time and gives the same expansion factor. However, prior to the fluctuation, the observers 
had different values of (p — (pi — 5 (p. 

It will also be convenient to allow a small interval dt for the time t\ at which the observers 
reach (p*. (For greater generality, I am using a time variable with an arbitrary parameter 
a). Then we are interested in the volume- weighted probability for an observer to reach 
(p = (pi — 5(p within the interval dt. This is given by \J(<fi — 8(p,ti)\dt, where J is the flux 
associated with the distribution V introduced in Section 5. The probability distribution for 
5(p is then proportional to 



<5p> /(Sp) 



rms 



6H*/H[ ~ (5p/p) rms < 1. 



(73) 




(74) 
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In the slow-roll range of ip, 



J « -(47r)- 1 ^ 1 ( V 9)^'(v?)Vi(v)exp(7it), 



(75) 



and we can use the approximate form of ipi(ip), Eq. (|36|) . Keeping only the factors dependent 
on Sip, we have 



dtP (Sip) oc exp 



2^ 



^-^(3-7i^r a ) 



d(Sip). 



(76) 



Note the a-dependence of the result. 

The average value of Sip can now be written as 



< Sip >- 



HI 
71 Hi 



3-d 



H„ 



(77) 



where the fractal dimension d is given by ( fi3|) and H max = max{H(ip)}. For a = 0, which 
corresponds to t = In a, 

(3 - d)H\ 



< Sip > 



a=0 



7lH[ 



(78) 



which is the same as our result (|7C|), apart from the numerical coefficient. But for a = 1, 
corresponding to the proper time t = r, the result is very different. The expansion rate Hi 
for the values of ip corresponding to the observable range of length scales is typically much 
smaller than the highest expansion rate H max , and thus 



<Sip > \ a= i 



TlH'i 



I 1 



(79) 



and 



< Sip > /(Sip), 



-_i — 2dH max Hi/ H[. 



(80) 



This is the result of Linde et.al. If one assumes that d ~ 1 and H max ~ 1, as it is often done 
in chaotic inflation scenario then it follows from Eq.fl^) that < Sip >^> (Sip) rms and 

thus < Sp >^> (Sp) rms . 
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VIII. DISCUSSION 



In this paper I have suggested a cutoff procedure which allows one to assign probabili- 
ties to different types of thermalized volumes in an eternally inflating universe. The main 
advantage of this procedure is that it is rather insensitive to the choice of the time coor- 
dinate t which is used to cut off the infinite volumes. In some cases there is essentially no 
t-dependence at all. Examples are the relative probability for the two minima for a potential 
shown in Fig. 4 (Section 6), and the probabilities for different fluctuation spectra (Section 
7). 

For a potential with multiple minima, we found that the main factor determining the 
relative probability of different minima is the slow- roll expansion factor (Section 4). This 
factor has a very sensitive dependence on the form of the potential V((p) in the corresponding 
slow-roll regions, and we expect the probabilities for different minima to be vastly different. 
If one finds that one minimum is more probable than another, this conclusion is not likely 
to be changed by a different choice of the time coordinate, even in cases where the results 
do show some t-dependence. Since this dependence is expected to be rather mild, it can 
affect only rare borderline cases where the two probabilities are nearly equal. 

One can accept this as a genuine uncertainty of the problem and make predictions only in 
cases where one minimum is much more probable than the others. An alternative attitude 
is to assert that there is a preferred choice of the time variable which gives the 'correct' 
probabilities. If this second approach is taken, then there is, arguably, a good reason to take 
a = (t = In a) as the preferred choice [|14" , |21 |. The only semiclassical variables that can 
be used as clocks in an inflating universe are the inflaton field <p and the scale factor a (at 
least in the simplest models). We are comparing volumes of the hypersurfaces <p = const, 
and thus the only remaining variable to be used as a clock is a. 

The cutoff prescription proposed in this paper can also be used to compare volumes in 
different, disconnected, eternally inflating universes. The world view suggested by quantum 
cosmology is that such universes may spontaneously nucleate out of nothing, and that the 
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constants of Nature may take different values in different universes ||22|| . The principle of 
mediocrity suggests that we think of ourselves as a civilization randomly picked out of an 
infinite number of civilizations inhabiting this metauniverse. The probability distribution 
for the constants may then be found by comparing the thermalized volumes in different 
universes. [Of course, one also has to compare the 'human factors' v civi see Eq. ([]])]. The 
volumes can be calculated along the same lines as in Section 4, except now the initial 
conditions are specified on the initial hypersurface at the 'moment of nucleation' and are 
determined by solving the Wheeler-DeWitt equation for the wave function of the universe. 

An important difference between comparing different minima of the potential and dif- 
ferent universes is that in the latter case the eigenvalues 71 (and 71) are not generally the 
same for the two universes. Since V*^ oc e - ^' ^ , the ratio Vi 1 ^ /Vi 2 ' ) will either vanish 
or diverge in the limit e — > 0, unless 71/71 = l[ 2 ^ /l^ ■ Hence, only the constants cor- 
responding to the largest possible value of 71/71 will have a non-zero probability. If the 
condition 71/71 = max does not determine all the constants of Nature, then the probabil- 
ity distribution for the remaining constants can be found by considering the e-independent 
factors in V* (such as [zi^] 3 ). These ideas were briefly outlined in Ref. |5[ and will be 
discussed in detail elsewhere. 
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FIGURES 



FIG. 1. 

Scalar field potential V(ip) with two minima and two slow-roll regions. The field values 
ip^ and y?! 2 -* indicate the end of inflation. In the region between ip^ and pif 1 the field 
dynamics is dominated by quantum fluctuations, and the slow-roll regions are bounded by 

ifQ^ and ip*\ 

FIG. 2. 

A slice through the spacetime of an eternally inflating universe. Thermalized regions of 
two different types are shown by different shades of grey. 

FIG. 3. 

A congruence of 'observers' ending up in different types of thermalized regions. The 
initial volume Vo is shown by a horizontal line, and the observer's wordlines by vertical 
lines. 

FIG. 4. 

A potential with a flat central regions bordering with two slow-roll regions. 
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